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SUMMARY 


There are three mechanisms by which a vertical magnetic field 
might be generated by a nuclear explosion. Asymmeiry of the nuclear 
device might produce a non-vertical dipole moment whose time variation 
would lead to an electromagnetic field with a vertical magnetic component. 
The magnitude of the vertical component due to this mechanism would 
depend sensitively on device design and orientation at the time of the 
detonation. 

Alternatively, Compton currents produced by explosion gamma rays 
might interact with the earth's field and generate a "back" field tending 
to cancel the earth's field in a region about the explosion. Such a 
canceliation would be equivalent to radiation from a Magnetic dipole and 
therefore could lead to a vertical magnetic field. 

Finally, on a much longer time scale, the blast wave motion of the 
heated (conducting) air surrounding the explosion might exclude the earth's 
field and thus generate a low frequency variation in the vertical component 
of the magnetic field. 

In this report only the latter two mechanisms are treated. It is 
Shown that due to air conductivity and the fact that gamma rays travel 
with light speed, the Compton currents do not exclude the earth's field. 
The blast wave motion does exclude the field and the magnitude of the 


resultant electromagnetic field is estimated. 
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THE ELECTROMAGNETIC SIGNAL DUE TO 
THE EXCLUSION OF THE EARTH'S MAGNETIC FIRLD 
BY NUCLEAR EXPLOSIONS 


There are two possible mechanisms by which the earth's magnetic 
field might be excluded from a large region about a nuclear explosion. 
First, the Compton-recoil electrons generated by explosion gamma rays 
interact with the earth's magnetic field, producing a "back" field. Approx— 
imate calculations described here indicate that this mechanism does not 
exclude the field. The other possibility is that the field is excluded due 
to the blast wave motion of heated (conducting) air. Calculations indicate 


that this mechanism is effective and its magnitude is crudely estimated. 


HOW THE GAMMA RAYS TEND TO EXCLUDE THE FIELD 

The "prompt" y rays emitted by a nuclear explosion generate 
Compton-recoil electrons in the surrounding air. The recoil electrons 
move roughly a few meters at sea level, predominantly in an outward 
radial direction. As the electrons move, they are turned by the earth's 
magnetic field, thereby producing an azimuthal current which tends to 


cancel the field in the region near the explosion, 


MECHANISM INHIBITING EXCLUSION 
ENT AL EE TUIN 

The slowing—down of the Compton electrons by the air produces 
ionization and makes the air conducting. The conductivity of the air 


leads to a conduction current which opposes the direct Compton current. 
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The question of whether the conduction current cancels the Compton current 
is complicated because the gamma pulse moves with the speed of light, 
which means that the displacement current cannot be ignored and Maxwell's 


equations must be treated in detail —- which is done in the next section. 


THEORY 
Denoting the conductivity by o and the Compton current by rs 


Maxwell's equations become 


> 109E . 44, > 
VxB= = a + ——(GE + j), (1) 
+ 4 aB 


g at adistance r (kilometers) from the burst point and atime t after the 


explosion is 


2 
e 

o(t, r) = or n(t, r), (3) 

c 

1 
where approximately 
—r /r t 1 
nit, r) %¥ 4x 107° —- ny [ at f(t! _rieje ht a (4) 
. xr Q 


In Eqs. (3) and (4), m is the electron mass, vo is the collision frequency 
of thermal electrons with air molecules, \ is the mean free path (kilometers) 


of the explosion gamma rays, Y is the explosion yield in kilotons, n is 


1 Multiple scattering of the gamma rays is neglected. 
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3 
D the fraction of Y appearing as gamma rays, Bf is the electron attachment rate 
to O75 , and f(t) is the temporal variation of the gamma-ray emission 
oO 
(J dt £(t) = 1). 
0 
The Compton current rf is 
. 9 elk t r—(t—t'v,. (tt )vo/h_, 
(t, r) x —2x10°e 5 ny f’ , at'e(t'— Ve vt—t'), (5) 
hr t— — 
Cc 
where R is the range of the Compton electrons and v is their velocity. 
Assuming that the electrons are scattered radially away from the burst 
point, 
> > 2t—-t'> 
x + 
V% vye. t+ Vg & (6) 


where the coordinate system has its origin at the explosion point with its 
polar axis along the earth's magnetic field and we consider only those 


. , 2 . 
electrons moving in the plane 60 = 7/2. v, is the speed of the Compton 


0 
mc v 


<= is the Larmor radius, B being the local 


electrons and p = y 
magnetic field. 

To obtain a precise solution of Eqs. (3) and (4) involves a prohibitive 
amount of calculation. Wecan, however, obtain a reasonably good 
approximate solution. For this purpose, we neglect effects arising from 
the spherical divergence of the Compton electrons and treat the local 


interaction of the electrons and the magnetic field as if the explosion were 


planar. This is justifiable as long as all spatial changes are confined 


“The "back'' field from electrons in other directions is less than that 
from electrons in the plane 9 = 7/2. 
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within a region whose radial dimension is small compared to the distance 
from the burst point. The results will be seen to justify this approximation. 


Within this approximation, 


— _ — + — 

E=E e, F e. (7) 

B-B) = b = be,, (8) 
where B, is the earth's magnetic field. Substituting Eqs. (7) and (8) into 


Eqs. (1) and (2) and using the planar approximation, 


ab 1 8E 4a | 

or o ant c (oF * Jy)» (9) 
aE 1 ab 

ya (10) 


From Eqs. (3) —(6), we find that o and Jy have the form 


o(t, r) = G(r) g(t—r/c), (11) 
J, {ts r) = G(r) h(t —r/c), (12) 


where G(r) varies slowly with r and may be treated as a constant. 


STEADY STATE SOLUTION. 
In view of Eqs. (11) and (12), b and E may be assumed to depend 


only on t = t~—r/e. Eqs. (9) and (10) then become 


_ ty = it) 4a[o(t) E(7) + 3,(7)], (13) 
dE(vr)  db(t) 
— ar ar 14) 


. 3The equation for F can be disregarded since it does not affect the 
B-field. 


RM-—2890 
D 
Solving Eqs. (13) and (14) leads straightforwardly to 
jy (7) 
E(t) = — “o(t). (15) 
and 
jy (7) 
b(t ) = Olt) + A, (16) 


where A is an arbitrary constant. Causality requires that for t < 0, 


E(t) and b(t) be zero. A Simple argument shows that 


| jy(™) 
am Olt) = 0, (17) 


Hence, A =Q and 


J 
b(t) = ? (18) 


Equations (15) and (18) for E(t) and b(t) satisfy Maxwell's equations 
and ali boundary conditions. They are therefore the required solutions. 

From Eqs. (5) and (6), it is seen that jy (7) < 0. Thus b(t) < 0 
and therefore adds to the earth's magnetic field. We have the surprising 
result that — not only do the Compton electrons not exclude the earth's 
magnetic field, but in fact they lead to an increase in the earth's field in 
the region where y(t) *# 0. This unintuitive result is a consequence in part of 
the fact that the electromagnetic field due to the displacement current 


moves with the gamma-ray pulse. If the current source had a velocity 


v << c, then the more familiar result that the field from a given current 
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opposes the field which generates it (Lenz's Law) would be satisfied (see 


Appendix I). 


Equation (18) is also a consequence of the assumed steady state for 


Jp 


current produces radiation which moves with the 


may therefore contribute to b. 


shown to be small. 


EFFECT OF TURNING ON j, 


In order to estimate the effect of turning on the current, we will 


assume oO is constant and J é 


and (10), 
ath 1 ab | 470 ab _ 4a 
ar? ce” at® 2 ot ¢ or 
Let 
. —ikr 
b = dre b(r, t) 
t., 
Then 
a 9b 2 2” A 
BB agg 2h 4 2c = — 4rick} 
2 ot b 
Ot 
and 
t , Af 2 2 D 
‘i ~ _ —_— — 
b = — 4rick f° dt e 21 o(t—7) sin (t-7)Vk c —(2r0)  & 
° Keo? — (200)" 


and o which ignores turning on the current. However, turning on the 


speed of light and which 


In the next section this contribution will be 


is turned on at t= 0. 


Then from Eqs. (9) 


o° 


(19) 


(20) 


(21) 


(22) 
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Hence, 
I 
t 00 0j,(r , 7) 
—27 oft — ; 
b(r, t) = ~4 f dte ( ” dr! 2 ({r—r'|, e(t—T)), (23) 
0) — 0 
where 
| — bir, 0) = 2B 2g 
and 


* 6) 


. 2 2 
i(|r—r'l, e(t—71)) =f dk cos k(r—r') sin e(t—T)V k -(Q@ra/c)" (24) 


0 4/ i? _ (Qno/c)” 
£3 (onio Ven ee) 0) ult le—r'Le). 


Alternatively, 


t 
9 _ 
b(r, t) = —2n far e Tot f(r tot-er, T) (r—ct+ecT, T) 
6) 


f r+c(t—t) 


| r. t 0 . 2 2, 2 
: _ d —~ V (t—71)°-—(r—r'! 
: Dan r jy , 7) a7 I,(2nio (t—T) —(r—r')}"/e )}. 
(25) 
Assuming again that 
jy(, t) = jp att —r/c) for t> 0, (26) 


co) 
where J dt g(t) = 1 and g(t) = 0 for t < 0, 
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we find 


j, g(t—r/c) t 
0 —2n ot — — 
b(r, t) = ————-———- (l--e ed ) — 27 j fat e an o(t—T} | g(2r—r/o~t) 
o 0 ‘ 
rt+c(t—T) 
| / 2 2, 2 
-f[ dr' g(t —r'/c) ~_ (2rio (t— T) —(r—r') ie ) | ‘ 
—- oe ae Or! O 
r—c(t—T) 
(27) 
The first term gives the steady state solution, Eq. (18). The remaining 
terms give the effect of turning on the current Jy: 
To illustrate the behavior of b we shall approximate g(t) by a delta 
function, 6(t). Thus 
j, &(t — r/c) _ _ 
b(r, t) = 0 -e anoty a er olt r/c) 
CO 0 
(ttr/c)/2 
_ —2ra(t—T) 9 ( alte wet Ne *) 
20 jo, dte a7 50 ania (t—t) —(r—ct) /e 
Or (28) 
j, 6(t—r/c) _ 
b(r, t) = v : (1 —e em Ot) 


(ttr/c)/2 _ _ 
— 27c jg — A dre 2m o(t—7) Jo (2ricy (tor)? =(r —e1) 7/07 ) 


0 


(29) 


It is straightforward to show that for o(t+r/c) >> 1 the last term of 


—27 ot 


Eq- (29) is Ofe ) (see Appendix II). Thus 
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b(r, t) = ——————_- + Of ). (30) 


More generally it can be shown that for an arbitrary positive and rapidly 
; —2n ot 
decaying pulse g(t), the last two terms of Eq. (27) are O(e ). 
Therefore, for constant 0, we have verified that the effect of turning 
on the current is small. 
On physical grounds this conclusion could have been expected since 
the radiation from turning on the current Je will be attenuated by 


propagating through the conducting region. We note, moreover, from 


Eqs. (23) and (24) and more directly from Eq. (29) that 
oa 
J ar blr, t) = 0 (31) 
— 0 


as expected from conservation of magnetic flux. 
Since no relevant new features are introduced by the spatial and 
temporal variation of the conductivity, corrections to Eq. (18) are 


expected in general to be small. 


ESTIMATE OF b(t) 


Substituting Eqs. (3) — (6) into Eq. (18) gives 


4 ; ; ; 
Provided, of course, that o is large in the region where the 
current is turned on. 
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10 
Vv 
I art e( + [rar tJ) elt Wo! (7 gy 
1 ‘co “t-Rie © 
O() = BT aye 7 
1 at 
2x10 [arte(r") eo BT T') 
0 
(32) 
To obtain a numerical estimate of b(t), we let f(1) ~ e, Yo ~~ ¢c, and 
yx 2. Then 
v ct/x 
+B 1- — 
b(n) = — (By + Ib(x)]) —2> S78 toe er) (33) 
4x 10 (£) l-e 
X 
a for 0< + < R/c and 
v R/x 
+ 1— — 
b(t) = — (By + |b(r) |) — = se) (34) 
4x10 (<) l—-e 
X 
: —8 , 
for Ric < tT. Ina time of the order 1/@ % 10° sec, b(t) increases 
from zero to its steady state value of 
“c a+B R/x 
b(t) = —(B. + |b(r)]) —S— 1-e'*(1—R/h) |. (35) 
Q 4 2 
4x10 (<) 
nS 
. il 8 
At sea level, v.~ 1.5x 10 per sec, B % 0.75 x10 per sec, and 
RX 3x 10° cm. Thus 
b(t) & — 0.033 (B, + |b(r)|) (36) 


or 


| b(t) ¥ — 0.033 Bp. (37) 


= 
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That is, from Eq. (8), 


B® —1.033 By ea- (38) 


HYDRODYNAMIC EXCLUSION OF THE EARTH'S FIELD 


Surrounding a nuclear explosion is an intensely ionized region in 
which the air conductivity is sufficiently high to freeze in the magnetic 
field. As a result, when the blast wave causes the air to expand, the 
magnetic field is carried along. At the same time the magnetic field 
which is external to the ionized region is unable to penetrate the ionized 
region and is pushed outward. 


The behavior of the magnetic field is described by 


i oB -vx(Yxp— V8 (39) 
c t Cc 4nro 

and 
vV-B= 0. (40) 


The details of the air heating and the hydrodynamic motion determine 
v and o as functions of time and space. These functions are so 
complicated that a precise solution of Eqs. (39) and (40) would involve a 
lengthy numerical calculation. Instead, we shall make an approximation 
which is adequate to provide semi—quantitative results. 

We assume that surrounding the explosion there exists an isothermal 


sphere of radius R(t) within which the air is highly fonized and the 
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relaxation time tT of the magnetic field is quite long, that is, 


rt aR) 


R(t) dRi(t) 
Rit) at ov 40 —— 


2 dt 
Cc 


>> 1. (41) 


Outside the isothermal sphere, the relaxation time is short and the 
conductivity can be neglected. With this simple model, the field at a 
distance r > R(t) from the explosion is 


2 3 
si yp > dl _ (RO) R(t) ine 
B = Bye, 5 By| B40) | ES sin 6 e, (42) 


2 3 
R(0) R(t) > 
~Bo|t- (Ro) | Lo COS 8 eps 


where again the polar axis coincides with the earth's field Bye, . 


Numerical results on the radiation and hydrodynamic flow for a 


sea level nuclear explosion give 


' 1/3 1/ 
R(t) & 5 + 120( 75) Y meters (43) 
1/3 
Y 
for 0 < ty il 3 < 0.05, where t is in seconds and Y in kilotons. For 
0.05 < ty! 
| 1/3 
R(t) = R ~ 50 Y meters. (44) 
max 


Equation (42) holds until Eq. (41) breaks down, after which the second 
and third terms of Eq. (42) decay to zero. The time constant for this 


decay 1s, approximately, 
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2 
Rnax 
410 —— , (45) 
c 
where o> is the conductivity of the isothermal sphere when R(t) = R nax . 
A rough estimate of a is 
Ty” 10+* sec . (46) 
2 
47 OOR max 2/3 
Consequently, a, is of the order of a few times Y seconds. 


Cc 


Since this relaxation time is so long, other phenomena, such as fireball 
rise and radiative and convective cooling, will determine the actual time 


for the magnetic field to relax. 
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APPENDIX I 
If in Eqs. (4) and (5), f = f(t—r/v), where v << c, then 


o(t, r) = F(r) g(t —r/v) ~ o(T) (I. 1) 


(t, r) = F(r) h(t — r/v) (tT), (I, 2) 


3 
" "6 
where tT = t—r/v. Substituting E = E(t) and b = b(t) into Eqs. (9) 


and (10), we find 


db v di anv 
aqE _ v db 
ayo oar (I. 4) 
Solving 
2 T 2 T 
E(r) = aaWvie) ip dt j (7 exp{ Sele) dt" as") (1.5)° 
1-(v/e)" “-0 1—(v/e)" “r! 
and 
b(t) = -— = E(r) + b(— 0), (I. 6) 


If b(—o) < 0 and Jy <0, then b(t) —b(—~) > 0; that is, the field 
behind the pulse is decreased in absolute value. Thus the field is 


partially excluded by the current — as expected from Lenz's Law. 


The non-vanishing of E(t) for large T is a consequence of the 
steady-state, one-dimensional character of the problem. 


ome 
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APPENDIX II 


Consider 


(ttr/c)/2 


dt eer alt—7) J (2nioV (t— 1) ~(r—et)“/e ) . (II. 1) 


—3 1 — Z). (II. 2) 


Then 


1 
_ 2 {etzle etoltale) p ac emotttrle)2 5 (anioV t2—(r/0)” yz). 
0 


(II. 3) 
For o(tt+tr/c)>>1, 


5) 


00 
I = ap Se TO fae e 7 Olttr/c)z I,(2nie Vt -(r/e)" fz ) 


+(ett+r) fe | (II. 4) 


But the integral in Eq. (II.4) can be evaluated exactly and leads to 


_ 9 1 —27 Ot 
I = afte + (ct + r) Ofe } . (IT. 5) 


Hence, 


_ 
) 
O 
_—~ 
av) 
~—_" 


(II. 6) 


